THEORY OF WAVEFORM-DIVERSE MOVING-TARGET SPOTLIGHT 
SYNTHETIC- APERTURE RADAR 



MARGARET CHENEY* AND BRETT BORDENt 

Abstract. Wc develop a theory for waveform-diverse moving-target synthetic-aperture radar, in the ease in whieh a single 
moving antenna is used for both transmitting and receiving. We assume that the targets (scattering objects) are moving linearly, 
but we allow an arbitrary, known flight path for the antenna and allow it to transmit a sequence of arbitrary, known waveforms. 

A formula for phase space (position and velocity) imaging is developed, and we provide a formula for the point-spread 
function of the corresponding imaging system. This point-spread function is expressed in terms of the ordinary radar ambiguity 
function. 

As an example, we show how the theory can be applied to the problem of estimating the errors that arise when target and 
antenna motion is neglected during the transit time of each pulse. 

1. Introduction. Synthetic-aperture radar (sar) images are formed using a moving antenna which 

transmits pulses of electromagnetic energy and measures the field scattered by objects (targets) in the envi- 
ronment. The measurements from different antenna locations are combined — typically by backprojection 
— and the resulting reflectivity map (target image) can display much finer resolution than that available 
from measurements at a single antenna location. (See, e.g. [1, 3, 4, 6, 8, 23, 25].) Standard SAR systems 
transmit the same high-range-resolution waveform repeatedly in order to measure the range (distance) from 
the antenna to various scene components. 

Below, wc develop the theory for a synthetic-aperture radar that makes use of a variety of different 
waveforms. We consider the case when there are multiple targets moving independently in the scene while 
the radar data are acquired. We assume that the targets are moving linearly, but we allow the antenna 
to fly an arbitrary, known, flight path and allow it to transmit a sequence of arbitrary, known waveforms. 
Our development is restricted to monostatic radar configurations, in which a single antenna is used for both 
transmitting and receiving. The problem is formulated in terms of forming an image in phase space, where 
the independent variables include not only the position of the target but also its vector velocity. 

A common approximation used in SAR imaging schemes assumes that the target is momentarily station- 
ary during the interval in which it interacts with an individual pulse (the so-called "fast-time" interval). In 
this "start-stop" approximation, the target is considered to move (relative to the radar platform) only during 
the "slow-time" interval between pulses. Because the speed of propagation of the pulse is much more rapid 
than the target motion, the start-stop approximation is justified for sufficiently short pulses. In our analysis, 
however, wc deliberately include the case of waveforms whose duration is sufficiently long that the targets 
and/or platform move appreciably while the data are collected. Figure 1.1 shows the approximate regime of 
validity of the start-stop approximation for for a system operating at 10 GHz (X-band). This figure plots 
V < A/T, where v is the relative velocity, A is the wavelength at the center frequency (3 cm), and T is the 
waveform duration. For shorter wavelengths, the curve moves towards the axes, and the region of validity 
is smaller. We see that the start-stop approximation is invalid for high-frequency systems, long-duration 
waveforms, or high- velocity targets. For example, a target moving 30 m/sec (67 mph) moves 3 mm in 100 
/xsec and 3 m in .1 sec. Such distances could easily be comparable to the radar wavelength, and the resulting 
errors would cause defocusing effects such as those seen in [7, 20] . This issue can become especially important 
when low-power, long-duration waveforms arc used. 

This paper is an extension of the work [2], in which we showed how to combine the temporal, spectral, 
and spatial attributes of radar data. In particular, the theory developed below shows how to combine fast- 
time Doppler and range measurements made from different spatial locations. This approach can be used, 
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Fig. 1.1. The approximate region of validity for the start-stop approximation for X-band (wavelength 3cm). The horizontal 
axis represents the waveform duration, and the vertical axis shows the relative velocity. For shorter wavelengths, the curve 
moves towards the axes, and the region of validity is smaller. 



for example, for SAR imaging when relative velocities are large enough so that target echo returns at each 
look are Doppler-shifted. Alternatively, this theory shows how to include spatial considerations into classical 
radar ambiguity theory. In addition, this approach provides a connection between SAR and Moving Target 
Indicator (mti) radar. 

SAR imaging of moving targets has been previously addressed using the start-stop approximation. In [7] , 
Fienup analyzed the phase perturbations caused by moving targets and showed how the motion affects the 
image. The work [18, 11] identifies ambiguities in four-dimensional phase space that result from attempting 
to image moving targets from a sensor moving along a straight flight path. The patent [10] and papers 
[12, 17, 20, 21, 30] all use the start-stop approximation to identify moving targets from radar data. The 
paper [9] uses a fluid model to impose conservation of mass on a distribution of moving scatterers, and a 
Kalman tracker to improve the image adaptively. 

The question arises of whether the Doppler shift can be exploited to provide extra information; the 
theory developed in this paper can be used to address this question. Some previous work [5, 14] suggests 
that Doppler information may be useful for obtaining velocity information about isolated moving targets. 
The present paper, on the other hand, handles an arbitrary distribution of moving targets as well as an 
arbitrary flight path. 

The validity of the start-stop approximation for SAR was previously considered [28] in the case of 
stationary targets, an antenna moving along a straight line, and chirp waveforms. This paper can be 
considered an extension of [28] to the case of moving targets and an arbitrary sequence of waveforms. 

In section 2 we develop a model for monostatic radar data that accommodates temporal, spectral, and 
spatial diversity. In section 3 we show how to process the signal to form a phase-space image, and we express 
the point-spread function of the imaging system in terms of the ordinary radar ambiguity function. In section 
4 we apply the theory to analyze errors due to the start-stop approximation. 

2. Model for radar data. The model for radar data involves a mmiber of ingredients: the model for 
wave propagation and scattering, the model for the incident field, the model for the received field, and a 
number of simplifying assumptions and approximations. We consider each in turn. 
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2.1. Scalar Wave Propagation. We model wave propagation witli the scalar equation 



(2.1) 



where c denotes the speed of light in vacuum and £{t, x) is one component of the electric field at time t and 
position X. Here s is a source term that consists of two parts: s = s™ + s^"^, where s'" models the source 
due to the transmitting antenna, and models the effects of target scattering. The solution £, which we 
write as 5*°*, therefore splits into two parts: 5*°* = S^'^ + S^''. The first term, satisfies the wave equation 
for the known, prescribed source s™. This part we call the incident field it is the field in the absence of 
scatterers. The second part of is due to the presence of scattering targets, and this part is called the 
scattered field. 

In scattering problems the source term s*"^ represents the target's response to an incident field. This 
part of the source function will generally depend on the geometric and material properties of the target 
and on the form and strengtli of the incident field. Consequently, s^'' can be quite complicated to describe 
analytically. Fortunately, for our purposes it is not necessary to provide a detailed analysis of the target's 
response and we note that for stationary objects consisting of linear materials, we can write s^'^ as the 
time-domain convolution 



where q{t,x) is called the reflectivity function. This function is conveniently described in the frequency 
domain. We write 



The frequency-domain reflectivity function Q(ijo,x) can display a sensitive dependence on oj. When 
the target is small in comparison with the wavelength of the incident field, for example, Q ^ uP' (this 
behavior is known as "Rayleigh scattering"). At higher frequencies (shorter wavelengths) the dependence 
on bj is typically less pronounced. In the so-called "optical region" (5(w,a;) is often approximated as being 
independent of a; (see, however, [22]); we use the optical approximation in this paper and simply drop the w 
dependence. In the time domain, this approximation corresponds to q{t,z) = S(t)Q(z). 

If the scatterer moves along the path x + F, where F passes through the origin at time t = 0, then the 
scatterer reflectivity translates as 



In general, both the source and the target will be moving; here we carry out all calculations in the rest 
frame of the stationary earth. 

2.2. A model for the field from a moving antenna. We model the antenna as an isotropic point 

source located at ^(t) (in the earth reference frame), transmitting (in its own reference frame) the waveform 
f{t). The corresponding equation for the time-domain field emanating from the antenna is then 




(2.2) 




(2.3) 



q{t,x)^6{t)Q{x-r{t)). 



(2.4) 



V^r" - Co-'a„£:™ = f{t)6ix - 7(t)), 



(2.5) 



where cq is the speed of light in vacuum. Consequently, the incident field is [14, 2] 



r'^{t,x) = - j 




t'-\x-y\/c) 
A'k\x — y\ 
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f{t')5{y-'t{t'))At'dy 



where we have used the Green's function for a point source [27] 

g(t,x) = ^ = — — dw, (2.7) 

' ' 4Tr\x\ J 8Tr^\x\ ' ^ ' 

which satisfies 

(V^ - co^du) g{t, x) = -S{t)d{x). (2.8) 
A more realistic antenna beam pattern can easily be included as in [3] . 

2.3. The Scattered Field. To obtain an equation for the scattered field, we solve (2.1) with the source 
on the right side given by (2.2). 

2.3.1. The Lippmann-Schwinger Integral Equation. We use (2.6), with / replaced by s^'', to 
express S^'' as a Lippmann-Schwinger integral equation 

£^'={t,x)= jj g[t-T,x- z)s^%T,z)dTdz 

= jj g{t-t',x- z)Q{z - T{T))£^°\t', z) At' dz. (2.9) 

2.3.2. The Born Approximation. For radar imaging, wc would like to determine Q from measure- 
ments of £^'^ at the antenna. However, it is difficult to do this directly from (2.9), because both Q and £^'^ 
in the neighborhood of the target are unknown, and in (2.9) these unkowns are mutliplied together. This 
nonlinearity makes it difficult to solve for Q. Consequently, almost all work on radar imaging involves making 
the Born approximation, which is also known as the weak- scattering or single-scattering approximation. The 
Born approximation replaces on the right side of (2.9) by f which is known. This results in a formula 
for £^'^ in terms of Q: 

f "^(t, x) « f|^(t, x):= jj g{t-t',x- z)Q{z - r(t'))^"(t', z) dt' dz. (2.10) 
In (2.10), we use (2.6) to obtain 

£%%t, x) = jg{t-t',x- y)Q{y - i{t'))£'%t' , y) dt' d\ (2.11) 

= - / '^^-''-\-y\i^^ f _ nt')) [ ^ft'-^" -i.-7(ni/c) . ,3^. 

J An\x-y\ J ^ "J An\y-j{t")\ ^ ' ^ 

2.4. The received field. Since we are modeling the antenna as an isotropic point source, the ffeld 
received at the antenna, which is located a,t x = j{t), is simply the field evaluated at the antenna, £g{t) := 

£h^{t,'rm 

We now make the change of variables y ^ z = y — T{t'), whose inverse is y = z -\-T{t'). This change 
of variables converts equation (2.11) into 

^"^'^ -'111 4.|7(i)-(z + r(iO)| ^^^^ A.\z + T(t')-^{t")\ " 

(2.13) 



Eq. (2.13) is a model for the field from a single localized moving target. The right side of (2.13) can 
be interpreted as follows. That part of the waveform / that is transmitted at time t" from location 'y{t") 
travels to the target, arriving at time t' . At time t' , the target that started at z is now at location z + T{t'). 
The wave scatters with relative strength Q{z), and then propagates to the receiver, arriving at time t. At 
time t, the receiver is at position 7(f). 

In the special case of a stationary reflector of unit strength positioned at the origin z = 0, we see that 
the times t" , t', and t are related by 

= t'-t"-\f{t")\/c 

= i-t'-|7(t)|/c. (2.14) 

If the scene involves a distribution of targets with reflectivities Qt, each moving on a different path F, 
then £g is integrated over all these target paths. To avoid dealing with such a path integral, we consider 
only the following simpler case. 

2.5. The Case of Pulsed Radar. For the rest of the paper, we assume that the waveform is a train 
of pulses and that the platform and target motion is slow relative to the speed of light. 

In particular, we take / to be a train of pulses of the form 

m = J2fm{t-T^), (2.15) 

where the delay between successive pulses is sufficiently large so that successive pulses do not overlap. 

A pulsed system allows us to introduce a "slow time" variable, namely the pulse number m. We reference 
the times t, t' , and t" to the following ("fast") times associated with the mth pulse. We write i/ = and 
t = for the solutions to (2.14) corresponding to t" = T^, so that 

= T:^-TZ-h{T^)\/c 

= r«-rl-|7(T^)|/c. (2.16) 

If the time intervals T^+i — are sufficiently short and the motion is smooth, we can make the expansions 

7(n=7(r,D+7(Tl)(t"-T,^) + --- 
7{t)=jiTi)+^iT^)it-T^) + --- 

r(i') = r(T4) + t{T;,)it' - t^) + • • • . (2.17) 

In making a distinction between the expansion times T^, and T^, we are including the case in which 
the platform is extremely distant from the target. A simpler version of the theory would be obtained by 
replacing the expansion times and by T^. 

2.6. Spotlight SAR. We now assume that the distance from the origin to the sensor position 7 is 

much larger than the distance from the origin to the target, and also much larger than the distance travelled 
by the target or sensor during any of the time intervals T^_|_i — T^. With these assumptions we can make 
the following expansions: 

\z + Tit') - 7(i)| = \z + r(T;) + t{T^){t' - T'J - 7(T«) - 7(T^)(i - T«) + • • • | 

= l7(T«)| - 7(T«) • [z + V{T'J + nrjit' - TU) - j{Ti){t - T^) 

\z + T{t') - 7(i")l = \z + nT'm) + nT'mW T'J ^(Tl) -i{Tl){t" - T^) + • • • | 

= |7(T^)|-7(T^^)- [^ + r(r;„) + r(T;,)(i'-r;„)-7(r^)(t"-T^) +••• (2.18) 
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Below we use the notation 7^ - 7(T^) , 7^ = 7(T^), r„ = r(r;„) , and = t{Tl^). 

To carry out the if' integration of (2.13), we use the approximations (2.18) together with the change of 
variables 



whose Jacobian is 



t" ^ ?' = t' - t" -\z + T{t') - 7(i")|/c 



at" 



dt" 



\dt"/dt"\ 



1 — . ' ^ / 



(2.19) 



(2.20) 



The t" integration of (2.13) contributes only when i" = 0, which occurs at a value of t" that we shall denote 
by T!^{t' ,z,Vm), which is given approximately by 



1 



1 1 ■ T I 



(1 + C • ^m/c) t' - YilMc + 7^ . (z + r„ - V^T'^ + 7^T^) /C 



(2.21) 



The time T^{t',z) is the retarded time, that is, the time when the field interacting with the target at time 
t' was transmitted. 

In (2.13) we also carry out the t' integration by making the change of variables 



whose Jacobian is 



t' ^i' = t-t' -\z + r(i') - 7(t)|/c 



dt' 



dt' 



\dt'/dt'\ 



-1 +7m • Vjn/C 



(2.22) 



(2.23) 



The delta function is supported only at i' = 0, which occurs at a value t' = T^{t, z, Vm) defined by 

= t-t' -\z + T{t')-'y{t)\/c. (2.24) 
With the approximations (2.18), we obtain 
1 



T^{t,Z,Vm) 



1 - 7m • Vm/C 

The expression (2.13) then becomes 



(1 - ■ ii/c) ( - bSl/c + ■ + r„ - o^n. + 7^32) /c] . (2.25) 



fm{(l>m{t,Z,Vjn)) Q{z) 



(4vr)2|7^ll7^l(l-7^-Wc)(l+/3D 



(2.26) 



where 



,{t,Z,Vm) ■■= Tl^{T^{t,Z,Vm),Z,Vm) - 
1 



([1 - & - R^iz, Vrr.)/c) - Rliz, v^) / c] - (2.27) 



and 



Pi 
Pi 



In 
^} 
In 



"Tm/C 



a,: 



1 + 7m • '"rn/c 

I /ml /m (.''^ m "m-^m^ Im-^mj 
l7ml ~ 7m ■ (-^ + ■'^m ~ '^'m^m 7m^m) * 



(2.28) 



The quantities 0^ and are determined by the squint angle (angle relative to broadside) of the antenna 
while it is transmitting and receiving, respectively, and a is the Doppler scale factor. 

We assume that the antenna illuminates a limited region of the earth and that all targets are contained 
in a bounded region in phase space. Moreover, we assume that the pulse repetition intervals are chosen so 
that the response from pulse m can be isolated from that of other pulses. Thus the data due to the mth 
pulse is just one term of (2.26): 



fm{<l>,n{t,Z,Vm)) Q{z) 



(2.29) 



Expression (2.29) applies to moving targets whose velocity during pulse w is t^m- To include multiple 
targets that may be moving with different velocities, we integrate over Vm- 



fm(i) 



/ Utt 



fm{(l)m{t,Z,V„i)) Q{z,Vm) 

(4^)^17^117^1(1- 7m •Wc)(l + /3^) 



■ d^^d^u,. 



(2.30) 



2.7. Example: Data model for circular SAR and constant-velocity targets. For the case of 
Circular SAR with no squint, we have 7^ • 7^ = = 7^ • 7^, which implies that (2.28) becomes 



/3m = = 

R^°{Z, Vm) = R-l^-{z + rm- VmT^) = R - • Z 
Rm°iz, Vm) = R - 7m • {Z + Tm - Wm^^) = R - • Z, 

where we have used the assumption that target velocity is constant to write 
(2.27) becomes 



C(i) = -J^^^ J fm{cl>°m{t,Z,Vm)) Q{z,Vm) d? Z <i?Vr, 

where from (2.27), the phase is 



(2.31) 

VmT'm- Thus (2.30) with 



(2.32) 



(j)°m{t,Z,Vm) = a^.,m {t - R^°iz,Vm)/c) - R1° {z , Vm) / C - T, 



1 + 



= t-Ti: 



R 



(?) 



R 

- -7m 



(f) 



.R B„ 



R z 

Z + Vm{t---'Jm-- 



(2.33) 



where Bm = 7™ + 7m the "bistatic bisector". For a system in which 7^ w 7^, Bm is simply the look 
direction. 



3. Image Formation. We process the data for each pulse, and then we combine the results over the 
pulse number m. 

3.1. Signal Processing at each Pulse. To process the data (2.30) at each pulse m, we apply a 
weighted matched filter, i.e., we correlate the scattered field (2.30) with a weighted version of the signal we 
expect to see from each point p and each velocity u. This results in 

Im{p,u) = - j {<l^m{t,P,Um)) (47r)2 1 1(1 ' ll ' t^m/c)(l + ^^)f^^(t) dt. (3.1) 

3.2. The Impulse-Response Function. The impulse-response function for each pulse is obtained by 
substituting (2.30) into (3.1). This results in 



-fm(P,Mm)= // Km(Um,P\Z,Vm)Q{z,V^)d.*VmA z 



where 



^m.{u^-,p--,Z,Vm) = y /m ('?!'m(i, P, Mm)) /m (<?i>m(i, 2:, ■«„)) 

In (3.3) we make the change of variables 
for which the Jacobian is 



1 _ / 

1 T' 
1 - 7m • '"m/C 



dt. 



dt 



dt' 



1 



I dt 



This converts (3.3) into 



Km{Um,P;Z,Vm) = Am { At(p, tt; Z, f) ) ^ 7m ^m/c 



1 - 7m • 1'm/c 



1 + 131 



where 



Am{<T,T) = J f:,{a[t-T])fm{t)dt 

is the wideband radar ambiguity function [26] and where 



ATto(p, Um;Z,Vm) 



^Rl,{p,u) Rl,{z 



ca 



+ (1 + /3^)T,^ 



v,m 

1 



a. 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



For the case of a stationary radar, (3.2) and (3.6) express the well-known fact [16] that a fixed radar 
can form a rangc-Doppler image of moving targets, and the system's ability to obtain range and velocity 
information is determined by the transmitted waveform via the radar ambiguity function. What is new 
here is that we have carefully accounted for the motion and spatial positions of the antenna and target 
distribution. 



3.3. Phase-Space Image Formation and Point-Spread Function. We now have an image at each 
puls{\ and a c;()rr(^sponding impulse-response function that depends on the pulse number m. This information 
could, for example, be fed into a tracking algorithm. In the present paper, we consider instead the formation 
of a phase-space image; in other words, we want to obtain information about both position and velocity of 
the targets. In order to combine information from multiple pulses, we make the assumption that the target 
velocity does not vary with m. We form the phase-space image as a weighted sum (3.2) over m, obtaining 

I{P,u) = ^/„(p,ti) 

m 

We write this as 

7(p,w) = J K{p,u;z,v)Q{z,v)d^zd^v, (3.10) 
where the full point-spread function K is the sum over m of the impulse- response functions (3.3): 

K{p, u; z, u) = ^ Km{p, u; z, v) 

m 

We write the point-spread function (3.11) as 



K{p,u;z,v) = ^ 

m 

Thus we see that the phase-space point-spread function (PSF) is a weighted, coherent sum of ordinary 
radar ambiguity functions, evaluated at arguments that depend on the difference in positions and velocities. 

This theory can potentially be applied to a variety of SAR image analysis problems involving moving 
targets, including the design of waveforms for tracking, analysis of resolution in phase space, etc. Below we 
show one example, namely an analysis of the validity of the start-stop approximation. 

4. Application: Validity of the Start-Stop Model. Almost all SAR signal processing and analysis 
makes use of the start-stop approximation, also called the stop- and- shoot approximation, which assumes that 
neither target nor antenna is moving while each interacts with the wave. 

4.1. The Start-Stop Model for Spotlight SAR Data. The usual start-stop model for monostatic 
SAR is 

where Rm denotes the distance from the stationary antenna to the target at z moving with velocity v. 

Below we use the theory developed above to extend the start-stop model to include the case of a fast- 
moving antenna very distant from the scene, in which case the antenna may not be at the same location 
when it transmits and receives. In this case, the start-stop approximation is to assume that the antenna 
is stationary while it is transmitting and while it is receiving, but it may move between transmission and 
reception. 

Under the start-stop assumption, (2.18) is instead 

\z + r(t') - 7(t)| « N + r{T'^) - 1{T^)\ « |7(T^)I - 7(T^0 • N + ^{T'J] + ■ ■ 
\z + T{t') - 7(t")l « k + nT'm) - l{Tl)\ « |7(r^)| - l{Tl) • [z + T{TU)] + ■■■. (4.2) 



E 



Km{p,u;z,v) 



Q?p,m(l f^m) 



Q{z,v)d^zd^v. (3.9) 



(1-/3^) 



(3.11) 



,ATmip,u;z,v) 



1 ^« / 

1 - 7m • U/C 
1 - 7m • Vjc 



(3.12) 



9 



Thus, using = VmT!^, the start-stop signal model is 



/ 



Q{z,V,n)d^Z&^Vm, (4.3) 



fm{t-{Rll''{z,v) + R^r{z,v))/c-Tl) 



(47r)2|7^ll7^l 
where 

Rll^^{z,v) = \^l\-^l.[z + vT'J 

R^^\z, v) = |7^| - 7™ • + ^T;;] . (4.4) 

Here we have used the same definitions for 7^,, 7^^, r/„, etc., as in sections 2.5 and 2.6. 
4.2. Start-Stop Image Formation for Spotlight SAR. The image is formed as 

7,,(p,«) = $^(47r)^|7^||7^| / fl{t-Tl-{Rli^\p,u) + R^^\{p,u))/c)£Ts{t)dt, (4.5) 

which gives rise to a PSF of the form 

K,,{p, u,z,v) = J2 f [t-Tl- (Rll^^ip, u) + R^'^iip, u))/c) 

U {t^Tl- {Rll'^z, v) + R^-\z, v))/c) dt 
= ^^„(1,AC), (4.6) 

m 

where 

Ar- = {Rli^^{z,v) + R^^\z,v))/c-{Rli^^{p,u)+R^^^{{p,u))lc (4.7) 

is the difference in travel times between the antenna and targets whose phase-space coordinates are (z, v) 
and {p,u). 

A focussed image will be obtained when 

= {Rli^^{z,v) + R^^^^{z,v)) - {R^^^{p,u) + R^^^^{{p,u)) 
= - 7™ • + ^T'J - (|7^|/c -ll-[p + uT'J) 

+ -ll-\z + vT'J, - (|7^|/c - 7™ • [P + uT'J) 

= ill + 7^) • [{P -z) + {u- v)T'J = Bm ■ [{p -z) + {u- vm , (4.8) 
where we have written 

Bm=7l+lt (4-9) 
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4.3. Mismatched Processing. If we incorrectly believe the start-stop data model to be accurate, we 
will form an image via the same processing as for the start-stop model, namely 



c 

where £^{t) is given by (2.30). Substituting (2.30) for £^{t) in (4.10), we obtain 

C'-(p,w) + i?^'-((p,«) 



f-(t)dt, (4.10) 



{P, U) 



f* i f -T^ 



f fm{4>m{t,Z,' 
J (1 - • V, 



The mismatched PSF is 



(l-7^.t;„/c)(l+/3^) 

Rlf\v, u) + i?«'^^((p, u) \ /^((/.^(t, z, -!;„)) 



(l-7^.^;™/c)(l + /3^) 
In (4.12), we make the change of variables t ^ t' = (pmit, z, Vm), whose inverse is 



(4.11) 



dt. (4.12) 



1 



it'+Tl) + 



c a. 



+ 



and whose Jacobian is 



dt 



dt' (l-/3«)a,,„,™' 

This change of variables converts (4.12) to 

1 + ^^ 



K^{p,u;z,v) / fm ( 



(1 - 



-(i'-Ar^) 



/m(i') 



(l + 7™-Wc)(l-/3^) 

1 



dt' 



(l + C-Wc)(l-/3^)' 



where 



At"" = T 



+ 



1 



(1 + /3^)C (! + » 



i?^'-(p.«)(l-/3«)a,^,„ i?^(z,t;) 



A focused image is the set of {p, u; z, v) that are at the ambiguity function peak for all m: 

= At^. 
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(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



The first line of (4.17) is the "Doppler" condition, and the second line is the "range" condition. We note 
that the reciprocal of the expression on the first line of (4.17) can be written 



(1 - P^)a^,m _ 1 - 7m • 7m/c 1 + 7m ' v/c 



(l - 7m • 7m/c - 7m * 7m/c) (1 + S„ ' V / c) 



(4.18) 



(4.19) 



= 1 - 7m • 7m/c - 7m ' 7m/c + • V / C + ■ ■ ■ , 

which, for the case of constant platform velocity (7 := 7^ = 7^), is 

^^^^^«l + Bm-(t;-7)/c+---. 

4.4. Error Analysis for Ideal High-Range-Resolution Waveforms. We expect that the start- 
stop approximation will be best in the case of a waveform whose ambiguity function is insensitive to its 

first argument, i.e., for a high-range- resolution waveform such as a very short pulse. For such an ideal 
high-range-resolution waveform, a focused image is attained when 



= Ar^ = T 



(1 - /3^)a„,. 



- 1 



+ (C'-(p,«)+ii^'-(p,«))(^ 



(1 + /3DC 



1 



(l + »- 

For the case of constant platform velocity, (4.20) becomes 

(1 l^rn)^v,m 



= T^ 

w rr 



1 



(l7ml - 7m • [P + wT^m] + l7ml " 7m " [P + "^1]) 



(1 - /3,f,)a„,, 



I lm\ im V ' im mj ' I /ml ! m V ' !m mj 



(T; - hl\/c)Brn ■ {v - 7)/c + (i?^^ -Brn-\p + uTj) (1 + B„ . (t; - 7)/c) /c 



(l -7m •7m/c) 



+ 



Rm Bm • Z 7m'7m(^m 17ml/'') 7m ' 7m (^m "I" 17ml/'') 

/c 



l7ml - 7m • + 7m(2;. + l7ml/c))J S„ ' t'/c 

= (TL - hl\/c)Bm ■ {v - 7)/c + {R^^ - B„ • [p + uTj) (1 + B„ • (t; - 7)/c) /c 

- (l - 7m • 7m/c) i?m^ - B„ • (z + 7T;) + 7^ • 7m(l7ml/c) - 7m " 7m(l7ml/c 
+ [l7ml - 7m • + "f^iTL + I7ml/C))] B„ • v/c /c 

= T'^B^ . {v - j)/c - |7^|B„ . {v - 7)/c^ 

+ {R^^ - -Bm • [P + «T;]) /c + (i?«^ - B„ . [p + wT;]) B„ . (t; - j)/c' 



R 



RT 



Bm-{z + jT^ 



/C + 7m • 7m/c 



b„.(z+7t;, 



(1 - 7m • 7m/c) 



7m • 7m(l7ml/c) - 7m * 7m(l7ml/c) 



-I- 



l7ml - 7m • + 7m (71 + |7ml/c))l B„ • v/c 



/C. 



(4.20) 



(4.21) 
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where we have written 

R^^ = hl\ + \i^\. 

We rearrange (4.21) as 

= Bm-[{z-p) + iv-u)T:^] /c 

- \jI\B^ . {v - 7)/c2 + «^ - B„, . [p + nr;]) B^.{v- j)/c' 



(4.22) 



R. 



RT 



Bm-{z+^T'^ 



7m • 7m(l7ml/c) " 7m * 7m(l7ml/c) 



+ [l7ml - 7m • + 7m (71 + |7ml/c))J Bm • v/c 

= B^.[{z-p) + {v-u)T'J Ic 

- |7^|Bm • {v - 7)/c2 + R^B^ • {v - ^)/^ - -Bm • [P + uT'J Bm-{v- 7)/c2 



~r 7m * 7m/'' 



i?: 



Bm-(;a; + 7?4) 



(l - 7m • 7m/c) 



7m ' 7m 



7m • 7m(l7ml/c) " 7m ' 7m(l7ml/c) 
|7^|Sm • v/c' - 7m • + imiTL + hm\/c)) Bm ' v/c 

hm\Bm • v/c - 7m • + 7m (^m + l7ml/c)) ^m " v/c 



(4.23) 



The first line of (4.23) is the condition (4.8) for focusing under the start-stop approximation. The bottom 
Une of (4.23) is of order {vjcf' and is neglected. Thus to leading order in {yjc), the error made in using the 
start-stop approximation is 

|7^|Bm • 7/c' - R^B^ • ^/c^ - Bm • [p + wT^] Bm-{v- 7)/c' 



"r 7m * '^ml ^ 



R^ -B^-{z^^T'J 



(l - 7m • 7m/c) 



7m * iml 7m * 7m/'' 



+ 



7^.(^ + 7^(T;, + |7S|/c))Bm-f/c^ 



(4.24) 



The third line of (4.24) is an overall image shift due to the sensor having movcid between the time of 
transmission and the time of reception. When this motion is negligible, we have 7^ — 7^ and B = 27, in 
which case the error (4.24) reduces to 



1 



4 (7m • [P + UT'^]) 7m • (t' - 7) - 27m ' 7r, 
+ 7m • (^ + im{T'm + \lm\/c)) 27m " ^ ■ 



(4.25) 
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If, in addition, there is no squint (7 • 7 = = i? • 7), the error (4.25) reduces to 



-B^ . [p + uTL] + 7™ • z\ ^^^^ _ ^/^^ ^ / -27,„-[p + ^7;g+7,„-^ \ (^^^ . ^/^^ (4_26) 



which vanishes for stationary targets. Thus for a stationary scene, we find that misfocusing due to the start- 
stop approximation can potentiaUy take place only in a squinted system. Below we compute the magnitude 
of this effect. 

4.4.1. EfTects on Image. Resolution in the time delay At is related to range resolution by At = 
2AR/c. Consequently, to compute the shift in the image due to the time-delay error (4.25), we multiply 
(4.25) by c/2. 

The image is affected only when the error due to the start-stop approximation is greater than a resolution 
cell, because in the backprojection process these contributions will not add correctly. It is well-known [3] 
that the SAR cross-range resolution S is determined by the aperture A9 and center frequency Ao as 

^ ^" (4.27) 



2A6' 



The angular aperture A9 is determined by the antenna path during which data is collected. Most systems 
are designed so that their down-range and cross-range resolutions are roughly equal. 

We consider the specific example of a low-earth-orbit satellite whose altitude is 300 km and whose 
velocity is 8 km/sec, traversing a straight flight path and forming a spotlight SAR image of a scene 500 km 
from the satellite. For such a rapid platform speed, the speed of most moving ground targets is insignificant 
in comparison and is consequently neglected in this example; however a similar analysis applies to the case 
of very rapidly moving targets. We use the assumption 7^ = 7^. 

We assume data collection starts at pulse m = and ends at pulse m = M. We denote the angle to the 
scene center relative to the antenna velocity vector by 9q for pulse m — 0, and 6m when for pulse m — M. 
See Figure 4.1. 




Fig. 4.1. Spotlight sar, geometry 



The duration of the data collection interval is denoted by T, which means for this case that the synthetic 
aperture is of length T\'y\, and the system transmits a pulse every T/{M — 1) seconds. To relate this pulse 
repetition interval, and hence the angular sampling interval, to the system's cross-range resolution, we apply 
the law of sines to the triangle shown in Fig. 4.1: 

sinAg sin 6*0 
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Solving (4.28) for T, we find that 



T = 



I 7m I sin Ag 
I7I sin^o 



25|7|sin6'o 




(4.29) 



where we have used the small-A^ approximation sin w A9, which from (4.27) can be written in terms of 
the system resolution as \q/{2S). 

In our simulations, we computed the error (4.25) for each pulse, plotted as a function of the resolution S 
and the starting angle ^o- In Figure 4.2, we show only the region in which the error is larger than a resolution 



Fig. 4.2. This shows the error as a function of squint and resolution for a target traveling lOOkm/hr in the direction 
opposite to the platform velocity (8 km/sec). We assume 3-cm wavelength (X-band). 

We see from Figure 4.2 that for a high-resolution system, when the squint angle ^0 is large, the start- 
stop approximation could lead to significant image distortion and smearing. This can be understood as a 
consequence of the difference in the way the Doppler scale factor a is treated in the start-stop approximation 
and the more complete model. Effectively, in the start-stop approximation, a is set equal to 1, whereas in 
the more complete model, we have expanded a in a series of powers of relative velocity. The error arises from 
the fact that we correlate an improperly scaled function f{at) against f{t). Contrary to common intuition, 
which is based on Doppler shifts rather than Doppler scale factors, the mismatch cannot be dismissed as 
a simple phase shift. In fact, however, this intuition argues that a very short pulse is sufficiently short 
that the start-stop approximation is justified, because the duration of its interaction with the target is 
"instantaneous". Such short pulses, however, are necessarily extremely wideband, and Doppler shifts are 
not appropriate descriptions of the Doppler effect. 

5. Conclusions. We have developed a theory that can be used to analyze waveform-diversity synthetic- 
aperture radar. This theory has many potential uses; we have shown, for example, that it can be used to 
analyze errors due to making the commonly-used start-stop approximation. The results show that even 
for very short pulses, these errors can be significant in squinted, rapidly-moving, high-resolution imaging 
systems. The errors are expected to be worse when longer, more realistic waveforms are used. 

The theory is potentially useful for analysis of imaging artifacts due to moving targets, and for tracking 
moving targets [9]. 



cell. 








squint angle (deg) 



error 5 (meters) 
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It is natural to hope that this theory can provide high-resolution phase-space images of multiple mov- 
ing targets. Wo note, however, that at each m, the radar ambiguity function of (3.12) provides us with 
two-dimensional information. The slow-time (position) variable m adds another dimension; thus we have 
three-dimensional data. In the case in which all positions and velocities are restricted to a two-dimensional 
plane, phase space is four-dimcnsional. Thus from three-dimensional data, we arc attempting to obtain 
four- dimensional information. This underdetermined imaging situation is currently poorly understood. In- 
vestigations of similar underdetermined imaging problems were carried out in [29, 19]; further exploration is 
left for the future. 
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